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Abstract 

Radiatively driven transfer flow perpendicular to a luminous disk was examined under a fully special 
relativistic treatment, taking into account radiation transfer. The flow was assumed to be vertical, and 
the gravity, the gas pressure, and the viscous heating were ignored. In order to construct the boundary 
condition at the flow top, the magic speed above the flat source was re-examined, and it was found that the 
magic speed above a moving source can exceed that above a static source 0.45 c). Then, the radiatively 
driven flow in a luminous disk was numerically solved, from the flow base (disk "inside"), where the flow 
speed is zero, to the flow top (disk "surface"), where the optical depth is zero. For a given optical depth 
and appropriate initial conditions at the flow base, where the flow starts, a loaded mass in the flow was 
obtained as an eigenvalue of the boundary condition at the flow top. Furthermore, a loaded mass and the 
flow final speed at the flow top were obtained as a function of the radiation pressure at the flow base; the 
flow final speed increases as the loaded mass decreases. Moreover, the fiow velocity and radiation fields 
along the flow were obtained as a function of the optical depth. Within the present treatment, the flow 
three velocity v is restricted to be within the range of u < c/\/3, which is the relativistic sound speed, due 
to the relativistic effect. 

Key words: accretion, accretion disks — astrophysical jets — radiative transfer — relativity — 
X-rays: stars 



1. Introduction 

Accretion disks are tremendous energy sources in the 
active universe (sec Kato et al. 1998 for a review). In 
particular, when the mass-accretion rate highly exceeds 
the critical rate, the disk local luminosity exceeds the 
Eddington one, and mass loss from the disk surface driven 
by radiation pressure takes place. 

Such a radiatively driven outflow from a luminous disk 
has been extensively studied in the context of models 
for astrophysical jets by many researchers (Bisnovatyi- 
Kogan, Bhnnikov 1977; Katz 1980; Icke 1980; Meha, 
Konigl 1989; Misra, Meha 1993; Tajima, Fukue 1996, 
1998; Watarai, Fukue 1999; Hirai, Fukue 2001; Fukue et 
al. 2001; Orihara, Fukue 2003), as on-axis jets (Icke 1989; 
Sikora et al. 1996; Renaud, Henri 1998; Luo, Protheroe 
1999; Fukue 2005a), as outflows confined by a gaseous 
torus (Lyndcn-BcU 1978; Davidson, McCray 1980; Sikora, 
Wilson 1981; Fukue 1982), or as jets confined by the outer 
flow or corona (Sol et al. 1989; Marcowith et al. 1995; 
Fukue 1999), and as numerical simulations (Eggum et al. 
1985, 1988). In almost all of these studies, however, the 
disk radiation fields were treated as external fields, and 
the radiation transfer was not solved. 

The radiation transfer in the disk, on the other hand, 
was investigated in relation to the structure of a static 
disk atmosphere and the spectral energy distribution from 
the disk surface (e.g., Meyer, Meyer-Hofmeister 1982; 
Cannizzo, Wheeler 1984; Shaviv, Wehrse 1986; Adam et 
al. 1988; Hubeny 1990; Mineshige, Wood 1990; Ross et 



al. 1992; Artemova et al. 1996; Hubeny, Hubeny 1997, 
1998; Hubeny et al. 2000, 2001; Davis et al. 2005; Hui 
et al. 2005). In these studies, however, the vertical move- 
ment and mass loss were not considered. Moreover, their 
treatments were restricted in the non-rclativistic regime, 
and the relativistic effects were not considered. 

In order to break such a situation, where radiation 
transfer has not been considered in the radiatively driven 
wind from the disk, we recently examined the radiatively 
driven vertical outflow - moving photosphere - in a lumi- 
nous flat disk (Fukue 2005b, c). In these papers, the radia- 
tive transfer flow was analytically or numerically solved, 
but the obtained flow speed was limited in the subrcla- 
tivistic regime, since the problem was treated up to the or- 
der of (v/c)^. In some astrophysical jet sources, however, 
the jet speed is mildly or highly relativistic. In gamma-ray 
bursts, the flow speed is supposed to be extremely rela- 
tivistic. Hence, in the next step, we should consider the 
transfer flow up to the full order of (v/c)'^ (cf. Fukue 1999; 
Hirai, Fukue 2001; Fukue et al. 2001; Orihara, Fukue 
2003; Fukue 2005a without transfer). Moreover, although 
it was incorporated in the previous paper (Fukue 2005b), 
which were up to the flrst order of the flow velocity, the 
effect of radiation drag must become much more impor- 
tant (Phinney 1987; Icke 1989). Other relativistic effects 
would further appear in the fully relativistic regime. 

In this paper, we thus examine the radiatively driven 
transfer outflow in a luminous flat disk within the fully 
relativistic regime of (y/c)^. At the preliminary stage, 
we ignore the gravity of the central object, the gas pres- 
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sure, and the viscous heating, and we treat the wind as a 
vertical one-dimensional flow without rotational motion. 

In the next section we describe basic equations in the 
vertical direction. In section 3 we examine the bound- 
ary condition at the flow top, and show the magic speed 
above the moving source. In section 4 we then solve and 
examine a radiative flow under the appropriate boundary 
conditions at the flow base and top. The final section is 
devoted to concluding remarks. 

2. Basic Equations 

Let us suppose a luminous flat disk, deep inside which 
gravitational or nuclear energy is released via viscous heat- 
ing or other processes. The radiation energy is trans- 
ported in the vertical direction, and the disk gas, itself, 
also moves in the vertical direction due to the action of ra- 
diation pressure (i.e., plane-parallel approximation). For 
the sake of simplicity, in the present paper, the radiation 
field is considered to be sufficiently intense that both the 
gravitational field of, e.g., the central object and the gas 
pressure can be ignored: tenuous cold normal plasmas in 
the super-Eddington disk, cold pair plasmas in the sub- 
Eddington disk, or dusty plasmas in the sub-Eddington 
disk. Internal heating is also ignored: the fiow in or near 
to the surface envelope of the disk. As for the order of 
the flow velocity v, we consider the fully special relativis- 
tic regime, where the terms are retained up to the second 
order of (v/c). 

Under these assumptions, the radiation hydrodynamic 
equations for steady vertical (z) flows are described as 
follows (Kato et al. 1998). 

The continuity equation is 



pcu = J const.). 



(1) 



where p is the proper gas density, u the vertical four 
velocity, J the mass-loss rate per unit area, and c the 
speed of light. The four velocity u is related to the three 
velocity v hy u = jv/c, where 7 is the Lorentz factor: 

7 = yrr^ = l/^y^-{v/c)^. 

The equation of motion is 

= H^^^ll^ [F^(l + 2«2) _ + Ph'u] , (2) 
az c 

where ^abs and Kgca. are the absorption and scattering 
opacities (gray), defined in the comoving frame, E the 
radiation energy density, F the radiative flux, and P the 
radiation pressure observed in the inertial frame. The first 
term in the brackets on the right-hand side of equation (2) 
means the radiatively driven force, which is modified to 
the order of u^, whereas the second term is the radiation 
drag force, which is also modified, but roughly propor- 
tional to the velocity. 

In the no-gas pressure approximation and without heat- 
ing, the energy equation is reduced to a radiative equilib- 
rium relation. 



0=j -CKabs-B7^-CKabs^'w^-|-2Kabs-f''7"7 (^^j du 

where j is the emissivity defined in the comoving frame. ut 



2Kabs-F7"7 



(3) 



In this equation (3), the third and fourth terms on the 
right-hand side appear in the relativistic regime. 

For radiation fields, the zeroth- moment equation be- 
comes 
dF 
dz 

(4) 



= P7 [j - CKabs-B + CKsce,[E + P)u^ 

The first-moment equation is 
dP p7 



dz 



-K,,,F{l + 2u^) + CK,,,{E + P)-fu]. (5) 
Finally, the closure relation in the inertial frame is 

cp(l + luA^cEa-luA+lF,u. 



(6) 



As a closure relation, the usual Eddington approximation 
in the comoving frame is adopted. Radiative quantities are 
then transformed from the comoving frame to the inertial 
frame, and we have the closure relation (6) in the inertial 
frame (see Kato et al. 1998 for details). 

Eliminating j with the help of equations (3), and using 
continuity equation (1), equations (2), (4), and (5) are 
rearranged as 

^du 

dz 
dF 

= (Kab. + Ksca)pu [c{E + P)-fU - F(l + 2^^)] , (g) 



CJ^ = (Kab. + A=sca)p- [F{1 + 2u^) - c{E + P)ju] ,(7) 



dz 

-( M \ ^ 

dz c 

Integrating the sum of equations (7) and (9) yields to 
the momentum flux conservation along the flow. 



= (Kab. + Ksca)p- [c{E + P)-fU - F {1 + 2u^)] . (9) 



cJu + P = K 



const. 



(10) 



In the subrelativistic regime, this relation is reduced to 
that derived in Fukue (2005b). Similarly, after some ma- 
nipulations, integrating the sum of equations (7) and (8) 
gives the energy flux conservation along the flow, 

cV7 + F = L (= const.). (11) 

In the subrelativistic regime, this relation means that the 
flux F is constant. 

At this stage, the basic equations are the equation of 
motion (7), the mass flux (1), the momentum flux (10), 
the energy flux (11), and the closure relation (6). 

Next, by introducing the optical depth r by 



dr = -(Kabs + Ksca)/Orfz, 

the equation of motion (7) is rewritten as 

cJ^ = -2 [F(l + 2u^) ~ c{E + P)ju\ 
dr c 



(12) 



(13) 



Furthermore, eliminating E with the help of equation (6), 
this equation (13) can be finally rewritten as 

7F(H-4u2)-4cP7M 



l-2u2 
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or 



F(l + 3/32)-4cP/3 



(15) 



1 - 3/32 

where (3 = v /c. 

We shall solve equations (14), (10), and (11) for appro- 
priate boundary conditions. Before this, we discuss the 
boundary conditions at the flow top in the next section. 

3. Magic Speed Above a Moving Photosphere 

When there is no motion in a luminous flat disk ( "static 
photosphere" ) , the radiation fields above the disk are eas- 
ily obtained. Namely, just above the disk with surface 
intensity Iq, the radiation energy density Es, the radia- 
tive flux i^s, and the radiation pressure Pg are (2/c)7r/o, 
ttIq, and (2/3c)7r/o, respectively, where subscript s de- 
notes the values at the disk surface. In the problem of 
radiation transfer in the accretion disk, these values are 
usually adopted as boundary conditions (e.g., Artemova 
et al. 1996). As will be shown below, however, the radia- 
tion fields above the luminous disk arc changed when the 
disk gas itself does move upward ( "moving photosphere" ) . 
Even in such a case, however, if the flow speed is small 
compared with the speed of light, the conditions for a 
static photosphere would be approximately adopted, and 
we can use these conditions in a previous paper (Fukue 
2005b), where the flow speed is limited in the subrcla- 
tivistic regime. When the flow speed is of the order of 
the speed of light, on the other hand, we should carefully 
treat the boundary condition for the moving photosphere. 
Thus, in the present paper, where the flow is treated in a 
fully relativistic manner, we must derive the exact bound- 
ary conditions above the moving photosphere. 

In addition, for the radiatively driven flow in the rela- 
tivistic regime, it becomes important the effect of radia- 
tion drag, which suppresses the jet speed (Phinney 1987; 
Icke 1989). For the flow above the static photosphere 
without gravity, Icke (1989) found that the magic speed oi 
jets becomes [(4 — •\/7)/3] 0.45 c. If the photophere is 
moving, however, such a magic speed will be also changed 
(cf. Fukue 2000). 

Hence, before we can examine the relativistic radiative 
flow, we must derive the radiation flclds above the moving 
photosphere and consider the boundary condition at the 
flow top (disk "surface"). 

Let us suppose the situation that a flat inflnite photo- 
sphere is not static, but moving upward with a speed Vs 
(= c/3s, and the corresponding Lorentz factor is 7s). Then, 
the direction and intensity of radiation are changed due 
to aberration and Doppler effects (cf. Kato et al. 1998; 
Fukue 2000). 

The transformations of the photon frequency v and pho- 
ton direction 6 between the inertial and comoving frames 
become 

^=7s(l-/3sCosg)= I — , (16) 

V 7s (H- /3s cos 6*0) 




0.2 0.4 0.6 0.8 1.0 



Fig. 1. Radiation fields above a moving photosphere with 
speed Va (= c/3s). Each component of the radiation fields is 
normalized by the values of a static case. Every component 
increases as the speed increases. 



cos9 



COS 6*0 -I- Ps 



(17) 

l-h^scos6io' ^ ' 

where subscript means the values measured in the co- 
moving frame and 7s = 1/ a/ 1 — /Jg • For incident radiation 
with ^0 = 7''/2 in the comoving frame, the direction cosine 
in the inertial frame is cosB = (3s. 

Furthermore, the transformation of the intensity / be- 
tween the inertial and comoving frames is 



/o = (^)*/=[74l-/3scos0)]^/ 



[7s(l + /3scos0o)]^ 



(18) 



Considering the Doppler effect (16) and aberration (17), 
the radiation fields above the moving photosphere in the 
inertial frame are calculated as follows: 

idn 





: 27r/o7; 



, 3 + 3/3s 



smt 



;i- /3s cos 6*)" 



= 2ttI, 



37s^ 



3 

- 37s Us 



0- 



/•cos ^ /3b 

Fs= / IcosOdn 
Jo 

2ttIq f™^"'* sin y cos 



(19) 



it Jo (l-^sCOSe)"* 



: 27r/o7; 



3 + 8(3s + 3/32 
6 



2ttIo 



372 +87sMs + 37/2 

6 



(20) 
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Fig. 2. Magic speed /3 above a moving photosphere as a 
function of the speed /3b of the photosphere (disk "surface"). 
When the photosphere is static, the magic speed is ~ 0.45c. 
When the photosphere is moving, however, the magic speed 
can exceed this limit. 



cP.= 



icos^edn 



2nIo f™" smOcos'^ Ode 



: 27r/o7,; 



2ttI, 



7s 



(1 — /3sCos^ 
1 + 3/3s + ■jpl 
3 

+ 37sMs + 



(21) 



That is, the radiation fields depend on the speed Vs of the 
photosphere, and every component increases as the speed 
increases (see figure 1). 

We must use these values of radiation fields as boundary 
conditions above the moving photosphere. 

In addition, the magic speed above the luminous infinite 
disk is obtained by the condition where the radiative force 
is balanced with the radiation drag force. Hence, from 
equation (2), wc have 

= F,(l + 2u2)_c(S3 + Ps)7w 

= [Fs(l + /32)-c(£;s + Ps)/3]7', (22) 

for the magic speed (3 (= u/7). 

Inserting equations (19)-(21) into this equation (22), we 
can obtain the magic speed above a moving photosphere as 
a function of the speed j3s of the photosphere. 



/3 = 



4-a/7 + 6/3s + (4 + %/7)4^ 
3 + 8/3s + 3/32 



(23) 



When the disk is static (/3s = 0) , this relation is reduced 
to that obtained by Icke (1989), /3 = (4 - ^/7)/3 - 0.45. 

In figure 2, the magic speed above a moving photosphere 
is shown as a function of the speed of the photosphere. 
When the luminous photosphere is not static, but moving, 
the magic speed can exceed the limit of ~ 0.45c, which Icke 
(1989) obtained for a static photosphere (cf. Fukue 2000). 



4. Relativistic Radiative Transfer Flow 

Now, we discuss our numerical solution of equations 
(14), (10), and (11) for appropriate boundary conditions. 

4-.1. Boundary Conditions and Singularity 

As for boundary conditions, wc impose the following 
cases. At the flow base (disk "inside") with an arbitrary 
optical depth tq (which relates to the disk surface density) , 
the flow velocity u is zero, the radiative flux is i^o (which 
is a measure of the strength of radiation field) , and the ra- 
diation pressure is Pq (which connects with the radiation 
pressure gradient and relates to the disk internal struc- 
ture), where subscript denotes the values at the flow 
base. At the flow top (disk "surface") where the optical 
depth is zero, the radiation fields should satisfy the values 
above a moving photosphere derived in the previous sec- 
tion. Namely, just above the disk with surface intensity 
/s, the radiation energy density Eg, the radiative flux Fs, 
and the radiation pressure Pg are, respectively, 



Ps = 2^Js 



cP, = 27rP 



37s^ 


+ 37s Us 


+ u2 




3 




37s^ 


+ 87sUs 


+ 3u2 
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(24) 
(25) 
(26) 



where Ug (= 7sWs/c) is the flow four velocity at the flow 
top and subscript s denotes the values at the flow top. 

Applying these boundary conditions to equations (10) 
and (11), we have two relations on the boundary values 
and mass-loss rate: 



Jc^Ws + cPs = cPo, 
Jc273 + P, = Jc^ + Po. 



(27) 
(28) 



Physically speaking, in the radiative flow starting from 
the flow base with an arbitrary optical depth tq , for initial 
values of Pq and Pq at the flow base, the final values of 
the radiation fields Ps, Ps, Ps, and the flow velocity Us at 
the flow top can be obtained by solving basic equations. 
Furthermore, the mass-loss rate J is determined as an 
eigenvalue so as to satisfy the bondary condition at the 
flow top (cf. Fukue 2005b in the subrelativistic regime). 

In the present full relativistic case, however, the final 
values of the radiation fields at the flow top depend on 
the flow velocity there, and the final values at the flow 
top cannot be analytically expressed by the initial values 
at the flow base. Hence, in this paper we determine the 
mass-loss rate as follows. 

In radiative flow with optical depth tq, we first give the 
final flow velocity Us (and 7s), instead of the initial value 
of Pq. Then, the final values of radiation fields Ps, Ps, and 
Ps can be fixed by equations (24)-(26). Next, we give a 
trial value for the mass-loss rate J, and the initial values of 
Pq and Pq can be fixed by equations (27) and (28). Since 
all the parameters are temporarily fixed, wc solve equation 
(14) from t = tq to t = 0. Generally, however, the obtained 
final velocity at r = is different from a given Ug. Thus, 
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we vary the value of J and follow iterative processes, so 
that the calculated final velocity coincides with a given 
final velocity Ug- 

Another point to be noticed is the singularity in equa- 
tion (14): i.e., the denominator of equation (14) vanishes 
when u = ±l/-\/2 or /3 = ±l/-\/3. This singularity orig- 
inates from the closure relation (6), and coincides with 
the sound speed of the relativistic (photon) gas. In other 
words, equation (14) has a form of the transonic wind 
equation for the relativistic (photon) gas with sound speed 
c/v^- In usual wind equations from a gravitating source, 
there exist transonic (critical) points, where both the nu- 
merator and denominator of wind equations vanish simul- 
taneously. In the present case, there also exist critical 
points, which yield that Uc = ±l/>/2 (/3c — il/VS) and 
Pc = sgn(?i)(\/3/2)Fc, where subscript c denotes the crit- 
ical point. In addition, from the linear analysis around a 
critical point, the velocity gradient near to the singularity 
is found to he du/dT\c= sgn(u)(3/4). 

We show an example of "critical solutions" in figure 
3, where the parameters are = 1, Fc/{ttIs) = 1, and 
J/ (tt/s/c^) = 1. As can be seen in figure 3, one of the crit- 
ical solutions (a solid curve) has a positive velocity, and 
is an outward breeze solution, which decelerates toward 
low optical depth. Another (a dashed curve) has a neg- 
ative velocity, and is an inward settling solution, which 
decelerates toward a high optical depth. Both solutions 
decelerate in the direction of the flow due to the radiation 
drag force. However, neither satisfies the present bound- 
ary conditions at the flow base and the flow top. Indeed, 
in order for flow to be accelerated, the velocity gradient 
du/dr should always be negative, but it is positive around 
a critical point for an outward solution. Hence, the flow 
in the present framework cannot pass through the critical 
point. Thus, for the present purpose, such "critical solu- 
tions" are inadequate, and the flow is always subsonic or 
supersonic in the sense that the flow speed is always less 
than or greater than c/a/3. And since we assume that the 
flow starts with u = at r = tq , we suppose a subsonic 
flow in the present paper: < l/V^ ~ 0.707. We shall 
discuss this singularity problem later. 

4-2. Subsonic Solutions 

Examples of the results for relativistic radiative flows 
in a luminous disk under the present boundary conditions 
are shown in figures 4 and 5. 

In figure 4 we show the final velocity Vs at the flow 
top (solid curves), the radiative flux Fq at the flow base 
(upper dashed ones), the radiative flux Fs at the flow top 
(lower dashed ones), the radiation pressure Ps at the flow 
top (dotted ones), and the mass- loss rate J (chain-dotted 
ones), as a function of Pq for several values of tq at the 
flow base. The quantities are normalized in units of c and 
Tr/g. For example, the unit of J is irls/c^. 

As can be seen in figure 4, as the radiative flux increases, 
the final flow velocity at the flow top increases, but the 
mass-loss rate decreases. Moreover, as can be seen in fig- 
ure 4, and similar to the subrelativistic case (Fukue 2005b, 
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Fig. 3. Example of "critical solutions". The parameters are 
Tc = 1, Fc/{itIs) = 1, and J/{-kIs/c^) = 1. These critical solu- 
tions pass through the critical points of the present problem, 
but do not satisfy the present boundary conditions. 

c), in order for flow to exist, the radiation pressure Pq at 
the flow base is restricted to be within some range. In the 
subrelativistic case without gravity and heating (Fukue 
2005b), the initial pressure Pq is proved to be restricted 
within the range of 2/3 < cPq /^s < 2/3 + to . In the present 
case, the initial pressure is also restricted to be within a 
similar range, but somewhat modified due to the relativis- 
tic effect: i.e., the radiative flux Fg is no longer constant, 
but depends on the flow final speed. At the one limit 
of Pq, the pressure gradient between the flow base and 
the top is maximum, the loaded mass diverges, and the 
flow final speed becomes zero. In the other limit of Pq, 
the pressure gradient vanishes, the loaded mass becomes 
zero, and the flow final speed becomes high. 

As already stated, in the present model, the flow final 
speed is supposed to be within the range of < Ug < 1/ 
(i.e., < Vs/c < l/-\/3)- As a result, using equations (25) 
and (26), the radiative flux Fg and the radiation pres- 
sure Pg are restricted to be within the following ranges: 
1 < Fg/inQ < 4.31 and 2/3 < cPg/{-Klg) < 3.73. Finally, 
although the flow final speed is less than the photon sound 
speed {vg/c < 1/a/3 ~ 0.577), the flow final speed can ex- 
ceed the magic speed above the luminous infinite disk 
K„agic/c=(4-V7)/3-0.451]. 

In figure 5 we show the flow three velocity v (solid 
curves), the radiative flux F (dashed curves), and the ra- 
diation pressure P (dotted curves) as a function of the 
optical depth r for several values of Ug at the flow top in 
a few cases of Tq. The quantities are normalized in units 
of c and nig. 

When the initial radiative flux Fq at the flow base is 
large, the flow is effectively accelerated, and the flow final 
speed becomes large. On the other hand, when the pres- 
sure gradient between the flow base and the flow top is 
large, the loaded mass J becomes large (cf. figure 4) . The 




Fig. 4. Final velocity Vg at the flow top (solid curves), the 
radiative flux _Fo at the flow base (upper dashed ones), the 
radiative flux Fb at the flow top (lower dashed ones), the 
radiation pressure Pg at the flow top (dotted ones), and the 
mass-loss rate J (chain-dotted ones), as a function of Pq for 
several values of tq at the flow base: (a) tq = 1, (b) tq = 10, 
and (c) TO = 100. The quantities are normalized in units of c 
and tt/s- That is, the unit of F and cP is wis and the unit of J 
is ttIs/c^. The left-side scale of the ordinates is for quantities 
except for J, while the right-side is for J. 



Fig. 5. Flow three velocity v (solid curves), radiative flux 
F (dashed curves), and radiation pressure P (dotted curves), 
as a function of the optical depth t for several values of Ub 
at the flow top in a few cases of tq. The parameters are (a) 
TO = 1, (b) TO = 10, and (c) to = 100, and from bottom to top 
of V, F, and P, the values of Us are 0.1, 0.3, 0.5, and 0.7. The 
quantities are normalized in units of c and Tr/g. That is, the 
unit of F and cP is it Is- 
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latter properties are somewhat complicated, and should 
be explained much more. 

When the optical depth tq is not so large (figures 4a and 
5a), the pressure gradient between the flow base and top 
is not so large. In this case, the pressure gradient becomes 
large, as the pressure Pg at the flow base becomes small. 
This is the reason that the loaded mass becomes large, 
since Pq is small (figure 4a). When the optical depth 
is large (figures 4b, 4c and 5b, 5c), on the other hand, 
the pressure gradient generally becomes large. Hence, the 
pressure gradient is large, as the pressure Pq at the flow 
base becomes large. This is the reason that the loaded 
mass becomes large, since Pq is large (figures 4b and 4c). 

Alternatively, in the less-luminous, sub/non-relativistic 
limit with small Ug, the radiation fields are roughly ex- 
pressed as 

F^F,, (29) 
cP ~ (2/3 + r), (30) 

which arc the usual Milne approximations (cf. Fukue 
2005b, c). Hence, in the less-luminous, small velocity case, 
the pressure gradient between the flow top and base be- 
comes large, as the optical depth is large. In the highly- 
luminous, relativistic limit with large Us, on the other 
hand, the flow velocity quickly reaches the equilibrium ve- 
locity, where the right-hand side of equation (14) or (15) 
vanishes, and becomes constant as 

where the radiative flux F and the radiation pressure P 
are also approximately constant according to equations 
(10) and (11). For example, in the limiting case of /3 = 
l/-\/3, cP/F = %/3/2. As a result, the pressure gradient is 
small even for large optical depth. 

In any case, the mass-loss rate J is not given arbitrarily, 
but is determined as an eigenvalue, similar to the subrel- 
ativistic case (Fukue 2005b). 

5. Concluding Remarks 

In this paper we have examined radiative flow in a lu- 
minous disk, while taking into account radiative transfer, 
in a fully relativistic manner (cf. Fukue 2005b, c for a 
subrelativistic regime). The vertical velocity u, the ra- 
diative flux F, and the radiation pressure P are numer- 
ically solved as a function of the optical depth r for the 
cases without gravity and heating. At the flow base (disk 
"inside") where the flow speed is zero, the initial optical 
depth, the initial radiative flux, and the initial radiation 
pressure are tq, Fq, and Pq, respectively; in the usual ac- 
cretion disk these quantities are determined in terms of 
the central mass, the mass-accretion rate, and the viscous 
process as a function of radius. At the flow top (disk 
"surface") where the optical depth r is zero, the radiation 
fields {Es, Fs, and Pg) should coincide with those above a 
moving photosphere with the final speed Vg with uniform 
intensity. In order to match this boundary condition, the 
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mass-loss rate J is determined as an eigenvalue. 

One of the relativistic manifestations is this boundary 
condition at the flow top. The radiation fields are not for 
the static flat source, but should be for the moving flat 
source. As a result, it is found that the magic speed above 
a moving photosphere can exceed that above a static pho- 
tosphere 0.45 c). 

Another relativistic manifestation is the existence of the 
singular point in the equation, which is related to the 
sound speed of the relativistic (photon) gas. Since the 
"critical solutions" for the present problem do not satisfy 
the present boundary conditions, we obtained subsonic 
solutions, which is always less than the relativistic sound 
speed of c/^/3. Nevertheless, we can find the relativistic 
solution, whose final speed is greater than the magic speed 
above a static source. 

The appearance of this singularity may come from the 
closure relation for the quantities of radiation fields, which 
assumes the usual Eddington approximation in the comov- 
ing frame (Kato et al. 1998). That is, in the comoving 
frame we assume 

Po' = ^^0, (32) 

where subscript denotes the quantities in the comoving 
frame. The appearance of the singularity suggests that 
the Eddington approximation would be violated in the 
relativistic flow, whose velocity is greater than c/^/S. This 
may be because the diffusion speed in the comoving frame 
cannot exceed c/>/3, and/or the diffusion would not be 
isotropic in the comoving frame. Indeed, if we assume that 
the diffusion in the comoving frame is no longer isotropic, 
but the factor 1/3 in equation (32) is variable, similar 
to the usual variable Eddington factor, we can formally 
extinguish the singularity. 

The plane-parallel assumption may also play a role in 
the existence of the singularity, as well as the existence of 
the radiation drag force, which is related to the closure 
approximation transformed to the inertial frame. The na- 
ture of critical points should be examined more carefully. 
In any case, the existence of the singularity does not mean 
that the flow velocity cannot exceed c/-\/3, but merely 
suggests a possible failure of the formalism adopted in the 
present analysis. 

The relativistic radiative flow investigated in the 
present paper must be quite fundamental to accretion disk 
physics, astrophysical jet formation, gamma-ray bursts, 
etc., although the present paper is only the flrst step and 
there are many simplifications at the present stage. 

At first, in a rigorous sense, the present situations, such 
as a plane-parallel assumption, are valid only in very spe- 
cific circumstances. Indeed, when the accretion rate ex- 
ceeds the critical one, the disk would puff up to be geomet- 
rically thick, and the plane-parallel approximation would 
be violated. In a rough sense, however, the present situa- 
tions could be approximately valid in such a thick disk, as 
long as the gas motion and the radiative flux are almost 
vertical. If the gas motion is not vertical, but expands like 
a spherical flow, the present treatment must be reexam- 
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ined. In contrast to the plane-parallel case considered in 
the present paper, the spherically expanding case is also of 
great interest to us, and should be examined in the future. 

In the present paper, wc have assumed that the opac- 
ity is frequency- independent (gray), and therefore the ra- 
diation pressure force is independent of the velocity. In 
a cold gas, however, atomic absorption opacities become 
important. Such opacities are frequency-dependent, and 
therefore the radiation pressure force would be a function 
of the Doppler shift between the main absorption troughs 
and the peak of the black body spectrum from the disk. 
Hence, the radiation energy flux L would generally be a 
function of the velocity u. In addition, in such a case of 
atomic opacities, the super-Eddington radiation-pressure 
force on the flow may arise from a larger absorption opac- 
ity. Since the acceleration nature would be drastically 
changed; such a line-driven radiative flow must be care- 
fully examined. 

Furthermore, we ignored the gravitational field pro- 
duced by the central object. This means that the flow 
considered in the present paper would correspond to nor- 
mal plasmas in the super-Eddington disk, pair plasmas in 
the sub-Eddington disk, or dusty plasmas in the luminous 
disk. In other cases, or even for normal plasmas in the 
super-Eddington disk, the gravitational field would affect 
the flow properties (cf. Fukue 2005c for a subrclativistic 
case). In particular, since the gravitational field in the ver- 
tical direction is somewhat complicated (e.g., Fukue 2002, 
2004), the infiuence of the gravitational field is important. 

We also ignored the gas pressure. This means that the 
radiation field is sufficiently intense. In general cases, 
where the gas pressure is considered, there usually ap- 
pear sonic points (e.g.. Fukue 2002, 2004), and the flow is 
accelerated from subsonic to supersonic. 

In the high-energy regime, including relativistic fiow 
and hot gas, the Compton heating and cooling would be 
important, although we have dropped such an effect. If 
the Compton heating of the gas by intense radiation fields 
works, the gas might be radiatively heated to very high 
temperatures, and the atomic absorption opacity could be 
significantly reduced. Such a radiative effect in the energy 
equation should also be examined in future work. 

Finally, it is left as an open question whether the sin- 
gularity appearing in the present paper can be passed 
through or whether the relativistic radiation transfer for- 
malism should be improved. 

The authors would like to thank Professors S. Kato 
and S. Inagaki for their useful comments and discus- 
sions. Valuable comments from an anonymous referee are 
also gratefully acknowledged. This work has been sup- 
ported in part by a Grant-in- Aid for Scientific Research 
(15540235 J.F.) of the Ministry of Education, Culture, 
Sports, Science and Technology. 
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